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STRESS- INTENSITY FACTORS FOR A WEDGE-LOADED 
EDGE CRACK IN A SEMI-INFINITE STRIP* 

by 

F. Erdogan and H. Terada** 

Lehigh University, Bethlehem, Pa. 


ABSTRACT 


The problem of a semi-infinite strip containing an edge crack is 
considered. It is assumed that the strip is loaded by a frictionless 
rigid wedge pressed into the crack. The resulting crack-contact problem 
Is formulated in terms of a system of singular integral equations. The 
behavior of the solution near the singular points is studied in detail. 

A series of numerical examples are given and the results are compared 
with those obtained by the method of boundary collocation and by the 
simple beam theory. 


1. Introduction 


Wedge loading of elastic materials is used in practice mainly in 
certain fracture toughness characterization tests and in fracturing 
solids by wedge-splitting or cleaving. In the former case the geometry 
of the medium is usually a rectangular strip which contains a syimetric 
edge crack. The specimen is known as the "compact tension specimen" if 
the length of the strip is of the same order of magnitude as its height, 
and as the "trousers' legs specimen" if its length is very long compared 
to itSheight. In fracturing of solids there is, of course, no standard 
geometry. However, in this case too a traction-free surface, an edge 
crack, and a bounded geometry with a rigid wedge pressed into the crack 
may be considered as the distinct features of the related solid mechan- 
ics problem. In solving these problems it is generally assumed that the 
external load is either a concentrated wedge force or is distributed in 
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some known fashion (see, for example, [1] and [2]). However, if the 
load is applied to the solid by pressing a rigid wedge into the crack, 
the problem is clearly a crack-contact problem in which the contact 
stresses and in some cases the contact area are additional unknowns. 

There are a number of solutions available in which the stress in- 
tensity factor is evaluated in the compact tension specimen and in the 
semi-infinite strip for various types of traction boundary conditions 
[ 1 - 4 ]. However, the related contact-crack problem which is described 
in Figure la does not seem to have been considered before. The contact 
problems considered in [5] and [6] discus: the wedging of infinite solids 
containing a V-notch prior to the formation of cracks. In this paper 
the wedge loading problem for a semi -infinite strip shown in Figure la 
is considered. The previous solutions for this geometry were based, 
in one form or another, on either a boundary collocation technique [1-3] 
or finite element methods (see, for example, [4] for a review). Here 
the elasticity problem is formulated in terms of a system of integral 
equations and the contact stresses and the stress intensity factor is 
obtained for various loading conditions. 

2. Formulation 


To formulate the problem first the infinite strip containing three 
symmetric cracks and shown in Figure lb is considered. It is seen that 
the formulation of the desired problem shown in Figure la may be obtained 
from that of Figure lb by letting d-Hi ans- c-K). In this study it is as- 
sumed that the wedge is rigid and the contact is frictionless. Thus, 
the problem is symmetric with respect to x and y axes in geometry as 
well as in loading. Consequently it is sufficient to consider one fourth 
of the region only (i.e., 0<x< e «and 0<y<h). In this region the x rnd y- 
components of the displacement vector satisfying the related equilibrium 
equations, the regularity conditions at x = <», and the boundary condition 
of zero shear traction on x = 0, 0<y<h (i.e., a xy (0,y) = 0) may be expressed 
as follows [7,8]: 
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.00 

u(x,y) * \ j ij{f ] + ^ 9] + syg 2 )sinh(sy) 

+ 1 (f 2 + Z+L + syg 1 )cosh(sy)}sin(sx)dx 


i r 4.1 - 1 s j x + isy 

+ ^ (^— ♦ 1 s | x ) g 2 ® ds » 

•oo 

s " I j { s 9 2 + SyQl )sinh(sy) 

+ j 9] + syg 2 )cosh(sy)}cos(sx)ds 


v(x,y) 


- 1 s | x + isy 
,e ds 


where < is the elastic constant (< = 3-4v for plane strain, i.e., for B»h, 
and (3-v)/(l+v) for generalized plane stress, i.e., for B«h, B being 
the thickness of the strip), and f 1 , f 2 , g ] , g 2 , and g 3 are unknown and 
are functions of s. These unknowns are determined by the remaining 
boundary conditions which may be expressed as 

a yy (x,h) = 0 , ° xy (x,h) = 0 , 0 <x<- , 

°xy(x,0) = 0 » 0<x< °° » 
o xx ( 0 »y) = 0 , c<y<d , 

> 

u(0,y) ■ 0 , 0<y<c , d<y<h , 
a (x,0) = 0 , 0<x<a Q , a<x<b , | 

JJ | 

v(x,0) = f(x) , a Q <x<a , 

— 7 v ( x »0) = 0 , b<x<* 



(3a, b) 

(4) 

(5) 
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Three homogeneous conditions (3) and (4) may be used to eliminate 
three of the unknown functions. The mixed boundary conditions (5) and 
(6) give a system of two dual integral equations to determine the remain- 
ing two functions. To obtain the relevant integral equations of the prob 
lem in this paper, however, a more direct method is used by defining the 
following two new unknown functions: 


jg u(0,y) = <t>(y) , c<y<d , 

(7) 

v(x,0) = *(x) , 0<x<b . 

(8) 

Note that outside the domains indicated in (7) and (8) the functions 4 >(y) 
and 4 >(x) are zero. Thus, differentiating (1), letting x = 0, and inverting 
g 3 may be expressed in terms of After eliminating three of the remain- 

ing four unknowns by using (3) and (4) and the stress-displacement rela- 
tions, the last unknown function may be expressed in terms t> and ^ by 
again differentiating and inverting equation (2). After some straight- 
forward manipulations the stress components involved in the mixed boundary 
conditions (5) and (6) may be expressed as 

,d 

* a xx (0 » y) = H 11 (y» r )<t>( r ) dr 

J c 


>b 

+ H 12 (y,t)i|/(t)dt , c<y<d , 

J o 

(9) 

* a yy(x*0) = j H 21 (x,r)<(>(r)dr 



fb 

+ [ H^ 2 (x,t)^(t)dt , 0<x<b , (10) 

■*0 


where y is the shear modulus, and the kernels H^j, (i,j = 1,2) are found 
to be 
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H n (y.r) ♦ r?y * h ll(y- r) - h ll l *’* r) 


(11a) 


, , C A,cosh(ys) + A,s1nh(ys) 

h n (y.r) * .1 . 2 — ds , 

1 J ° 1 ♦ 4she- Zsh -e- 4sh 


(lib) 


A^ (s) * 2e‘^ 2h " r * s (-2s 2 h 2 + 2s 2 rh - 3sr 


- sre 


" 2sh + 4sh - 2 - 2e' 2sh ) , 


(11c) 


A 2 (s) = -2sye- (2h - r)s (l - 2sh ♦ 2sr + e - 2sh ) , (lid) 


u t +\ - 2t (t 2 -y 2 ) 9 

(y z +t* ) 2 2 


® B^coshCys) + B 2 sinh(ys) 

-4sh 

o 1 + 4she - e 


sin(tsjds , 
(12a) 


B-j ( s ) - e’ 2sh (l + 2hs - 2h 2 s 2 - e' c,ls ) , 


-2hs< 


(12b) 


B 2 (s) - e’ 2hs ys(l + 2hs - e"“ ,:> ) , 


-2hs< 


(12c) 


, , 2r(r 2 -x 2 ) . „ f" C cos (sx) ds 

H 21 ,X * r) = * J 0 ! + 4hse- 2hs e- 4hs 


(13a) 


C( s ) = e" 2hs [e sr (2h 2 s 2 - 2hs 2 r - sr + sre” 2hs ) 


- e 


-sr 


-2hs- 


(2h 2 s 2 + 2hs 2 r + sr - sre‘ )] , 


(13b) 


72 


(«.t) ■ F? + tk + h 22 (x>t) ’ 


h 22 (x»t) = 2 


» _-2hs^ e -2hs 
o 1 + 4hse' fc,,J - e 


(14a) 


r e ~ 2 [) 2s2 - ~ 2 jj s ~ ^ [sin s(t-x) 

-2hs -4hs L 
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+ sin s(t+xjds . 


(14b) 


For c>0 and d<h it can be shown that the kernels h^ , H-| 2 , H 21 » and 
h 22 are bounded within the rectangular closed intervals in which they are 
defined. In this case, referring to Figure lb, the definitions (7) and 
( 8 ), and the mixed boundary conditions (5) and ( 6 ), it may be seen that 
the unknown functions <t> and 4 / must satisfy the following conditions of 
single-val uedness of displacements: 


»d 

<?>(y)dy = 0 
J c 

’b 

iKx)dx = 0 
-b 


(15) 

(16) 


Even though this problem is not very practical, it can be solved without 
any difficulty. For example, if in addition to wedging, the strip is sub- 
jected to a uniform stress = a Q , the integral equations (9) and 

( 10 ) may be expressed as 


f H n (y,r)$(r)dr + H, ? (y,t)ih ( t)dt + ^ 
Jc Jo J a 


1+K 

= - n - 3 — a - 1 

4u 0 J 


H 12 (y,t)f(t)dt 


H 21 (x,r)<f>(r)dr + f 0 H 22 (x,t)t |> 1 (t)dt + j" 
c J 0 ^ 


f d 

H 22 (x,t)f(t)dt , 0<x<a Q , 

J ao 


f H 91 (x,r)<t>(r)dr + [ H 99 (x,t)>j » 1 (t)dt + j H 

J c Jo Ja 


H 22 (x,t)f( t)dt , a<x<b , 


where it is assumed that 


12 (y»tH 2 (t)dt 
c<y<d , (17) 

H 22 (x,t)v 2 (t)dt 

(13) 

22 (x,t)v 2 ( t)dt 

(19) 
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< 


( 20 ) 


v(x,0) = *(x) 


4 , i ( x) , 0<x<a Q 
f(x) , a Q <x<a 
( x ) , a<x<b 


and f ( x ) is known from the wedge profile. The unknown functions $, 4 ^ , 
and 4*2 are obtained from the system of singular integral equations (17-19) 
in which only the main diagonal terms in the kernel matrix contain Cauchy 
singularities. For the integral equations defined on a union of arcs this 
scheme of increasing the number of the- unknown functions and integral equa- 
tions in such a way that each equation is defined on a single arc makes 
it possible to solve the system numerically in a straightforward fashion. 
Otherwise, due to the complicated nature of the related fundamental func- 
tions, it would have been difficult and quite cumbersome to obtain a nu- 
merical solution. The index of (17) is +1 and the conditon (15) is still 
valid. The indexes of (18) and (19) and the corresponding additional 
conditions depend on the wedge profile. 

3. The End Problem 


If we let d=h and c = 0 the problem becomes a stress-free "end prob- 
lem". In this case it can be shown that for y = h * r the kernel h-j] given 
by (lib) is not bounded. To show this and to separate the singular part 
of this kernel consider the leading exponential terms in A-j(s) and A 2 (s), 
which may be written as 

A, (s) = 2e~ (2h ” r)s (-2s 2 h 2 + 2s 2 rh - 3sr + 4sh - 2) , 

I 00 * 9 

A 2 Js) = -2sye' (2h " r)s (l - 2sh + 2sr) . (21a, b) 

By adding and subtracting A^ and und e r The integral in (lib) and by 
using the relation [10] 
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f,ai.(2h-r)s <f s1nh(sy)1 > ds = r e -(2h-r)s fslnhtsy)^ ds 
j s e (cosh(sy)J dr ^ Jo (cosh(sy) J 


_d^ 

dr 111 


i J y I . 
P RI ^VrjJ ’ 


( 22 ) 


the kernel h^-j(y,r) may easily be expressed as follows: 


h n (y,r) = h $ (y,r) + h b (y,r) , 

h (V r ) « 1 + 6111 - 7 ) + /( 6 -y)^ 

V y,r) 2h-y-r (2h-y-r)^ (2h-r-y) J * 


(23) 


“ [-1 + 6(h-y) £ - 2(h-y) 2 1 ^y , 

f oo A ( 3 ) 

o t[ H4she- l2sh ' -- 4sh ‘ A l» (s)]c0sh(ys) 


(24) 


-e 


+ [ 


A 2 (s) 

^rr ^rr - A 2 Js)]sinh(ys)}ds 

1+ 4she ‘ sn -e ^ sn L 


(25) 


The kernel h b is bounded in o<(y,r)<h and the singular part h s together 
with the singular terms l/(r-y) and l/(r+y) in (11a) constitute a general- 
ized Cauchy kernel. 

To investigate the behavior of the solution of (9) and (10) at the end 
points (x = 0, y = h) and (x = 0, y = 0) let 


4>(r) = 


Hr) * f(r)e ,1B 
r“(h-r) 6 r°(r-h) e 


4>(t) = 


9(t) , 9(t)e m 

t 3 (b-t) T t°(t-b) Y 


0<r<h 


0<t<b 


(26) 

(27) 


where f is Holder-continuous in 0<r<h, r a (r-h) s is a definite branch vary- 
ing continuously in 0<r<h, and it is assumed that 
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f 


I 


1 


! 


I 


f(0) t 0 , f(h) f 0 , 0< Re(a,e) < 1 


(28) 


Similar statements hold for g(t). At the common end point (r=0, t = 0) 
the functions $ and i|> must have the same singular behavior, hence the 
common power a in (26) and (27). Define now the following sectional ly 
holomorphic functions: 


rh 


F < z l>=7 


ilnl 

r * z i 


dr 



(29a, b) 


where the variables y and r are on the line Re(z-|) and x and t are on 
Re(z 2 ). Analyzing the asymptotic behavior of, for example, F-j (z) around 
the end points z-j = 0 and z^ = h by following the procedure described in 
[11], we obtain 


Hz } ) = 


f(0)e 1,ia 1 . ft hj ] +F(Zl) 

h 8 sinna Z 1 h a sin*B ( z-j -h ) e 0 


(30) 


h iirl dr = f (0)cot»a J_ . f(h)cotnB 1 + c ( y ) 

* o r -y h 8 y a h a (h-y ) 6 1 




I [ h t)r)dr . F(2h ) fill! ] — - ♦ F,(y) 

* Jo r -(2h-y) (.“sinus {h . y) 6 3 


, (31) 


(32) 

(33) 


where the functions Fj , (j = 0,..,3) are bounded everywhere witn possible 
exceptions of the ends o and h where they may have singularities weaker 
than zj a and ( z-j -h) ” e , respectively. The remaining two terms coming 
from (24) may be obtained from (33) by differentiation. Similarly 

g(z 2 ) - aMel“_L. _aM +g 0 (z 2 ) , (34) 

b Y sinna z 2 b slnity (z 2 -b) Y 
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I f b till dt . g(0)cotwa _1_ _ g(b)cotTiY _J + G , s (35) 

* J 0 t-X « b Y x“ b“ ( b . x) Y V X ' ’ W 


l 

TT 


’b 

o 


tin 

t+x 


dt 


G(-x) 


9(0) 

b Y sinna 


^ +G 2 (X) 


(36) 


where G c , G-j , and G 2 are again either bounded or have weaker end singular- 
ities. 

Consider now the singular kernels in and Hg-j given in (12a) and 
(13a). Note that these kernels are identical and may be expressed as 


2 + 2 y ( fV ) 


(37) 


Now, observing that 


2t _ 1 1 

t z +y 2 " t-iy t+iy * 


and at Z 2 = + iy G is holcmorphic, we obtain 


^ jzhrz ^(t)dt = G(iy) + G(-iy) 
0 L y 


2g(0)cos_(ja_/_2]_ 1 + G (y) f 
b Y sinna t 


(38) 


£ 1 f XFW * (t)dt * - 23(P ! COS(TO/ - 2) 5 ^ My) ,(39) 


dy 7T J 


. Y • 4 

b sinira j 


where, again at y=0 G 3 and G 4 are either bounded or have singularitites 
weaker than y‘ a . Similarly 


2 

IT 


■h 

0 


r(r 2 -x 2 ) 
Cr^T Jt 


4 »(r)dr 


2f(0)cos(W2) 1-g 
h^sinna * a 


(40) 
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Note that in 0<y<h a xx (0,y) is bounded and with 4 and ip as de- 
fined in (26) and (27), the integrals involving the Fredholm kernels in 
(9) and (10) contribute only bounded terms. Thus, substituting from (31), 
(32), (33), (24), (38), and (39) into (9), multiplying both sides first 
by (h-y) 6 and letting y-»Ti, and then by y a and letting y-*0, we obtain the 
following relations: 

[2(1-B) 2 - (1+COSnB)] • 0 , (41) 

h sinus 

iiOi (cot™ ♦ -J— ) + alfil 2(l-o) . 0 . (42) 

Sinua t)i Sinna ' 

The characteristic equation (41) giving B is identical to that of edge 
crack and 90-degree infinite elastic wedge [8]. Note that (41) ha", no 
root in the strip 0<Re( & ) <1 - Hence the function <$>(r) has no power sing- 
ularity at r = h. 

Similarly, substituting from (40), (35), and (36) into (10) we obtain 


9^1 cotiTy = 0 , 

b a 


(43) 




cos ( rcct/2 ) + g(0) 
S i nira (jY 


(cotua + 


_J 

Si nira 


) = 0 


(44) 


Since g(b) + 0, (43) gives y = 1/2 which is the standard value for an 
imbedded crack tip. Also since f ( 0) and g(0) are finite and unknown, 
(42) and (44) gives 


1+COSTTa 
sin ' 2 -it a 


[(1 + COSua) - 2(l-a) 2 ] = 


0 


(45) 


As expected this is essentially the same characteristic equation obtained 
from (41) for the corner (x = 0, y = h). Hence, at the corner (x = 0, y = 0) 
the functions <t> (y ) and ^( x ) have no power singularity. 


We can now proceed and show that <f>(r) is a bounded function in the 
closed interval 0<r<h, and ^(t) has the following form 

*(t) = , 0<t<b , (46) 

(b-t) Y (t-b) Y 

where g is bounded in 0<t<b. Starting with the definitions (29) the 
following asymptotic relations may be obtained [11]: 


F( z-j ) = - 1 4>(0)logz 1 + | 4 >(h)log(z 1 -h) + (^(z^ , 


(47) 


r h 


fh 


iill 

r-y 


dr = - 1 4>(0)logy + 1 *(h)log(y-h) + C^y) , (48) 


(49) 


o dr " ‘ 1 4>(0)l°g(-y) + C 2 (y) , 

G(Z2 ) = .iaM logZ2 .i^_l_ + c 3 (22 ) , (go) 

M b W dt = -^ l09x - 9(b)c0t ^ 7 T^V +c 4 (x ) 


(b-x) 


; vp ( t ) dt = G(iy) + G(-iy) = - | logy + C 5 (y) , 

b (52) 

y j o tjgr *U)dt • c 6 (y) , 

~ f jTZ+yz 4 > ( r ) dr = - — <j>(0)logx + C^(x) » 

* n 


1 f 4^- dt = - 1 3(01 i 0 g(. x ) + c 8 (x) , 


TT ]q t+X 


* b 


(53) 

(54) 

(55) 


where C ,..,C n are bounded functions, 
o o 

In equations (9) and CO) the left hand sides are bounded and the 
terms involving Fredholm kernels give finite quantities. Thus, substi- 
tuting from (47 ) - ( 55 ) into (9) and (10), combining the finite quantities, 
and observing that log(-z) = logz + iri , we obtain 
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* 


I 


F 


- | C4>(0) + ^-j^-Jlogy + \ 4>(h)[log(y-h) - log(h-y)] 


+ C 9 (y) = C 1Q (y) 


(56) 


- 7 [*(0) + ^f^logx - g(b)cotiTY 

7T b T 


(b-x) Y 


+ C-j-j(x) = C-j 2 ( x ) , 

(57) 


where C 9 ,..,C 12 are bounded functions. Since g(b) is finite (57) gives 
cotuy = 0, and y = l/2, which is the standard result. The logarithmic terms 
in the second term of (56) cancel. For the remaining unbounded terms in- 
volving logy and logx to vanish we must impose the following condition: 

♦ (0) + ^ = 0 . (58) 

b Y 

Of course, (58) is nothing but the statement of the physical condition 
that, in the absence of shear tractions, after deformations there is no 
change in the -90-degree angle at (x = 0, y = 0). This follows from 

♦ (y) ■ u(+0.y) , *(x) - £ v(x,+0) = SteL. . (59) 

3X (b-x) 

In solving the integral equations for the free corner problem with <j> and 4/ 
as the unknown functions (i.e., for a Q >0, c = 0 in Figure 1), (58) must be 
used as an additional condition. 

Referring to Figure 1 if a Q = 0, around the end point x = 0, y = 0, the 
function 4 /(x) is known and, since the wedge if frictionless, the contact 
stress is a bounded normal traction. Therefore, the analyses given in 
this section, including the condition (58) regarding the behavior of 4>(y) 
remain to be valid. 


4. Exampl es 

As a first example consider the (frictionless) wedging problem shown 
in the insert of Figure 2 where it is assumed that c s 0, d = h, and the 

-’ 3- 


I 


I 


I 


I 


wedge is "flat", i.e.. 


— v(x,0) = f(x) = 0 , 0<x<a . (60) 

Referring to Figure lb, it is also assumed that a Q = 0. This assumption 
is, of course, subject to the verification that in the region (0<x<a, y=0) 
the contact stress ayy(x,0) is compressive. In this problem the tractions 
on the crack surface are zero, the external load is provided by the wedge 
which has a total thickness 2v Q , and the integral equations (9) and (10) 
are still valid and may be expressed in the form (see equations (9) and (10)) 

fh fb 

j H-j 1 (y,r)4»(r)dr + j H 12 (y,t)^(t)dt = 0 , 0<y<h , (61) 

rh rb 

j H 2 i (x,r)<f>(r)dr + J H 22 (x,t)^(t)dt = 0 , a<x<b . (62) 

The unknown function <t> is bounded in 0<y<h and is zero at y s 0, whereas \p 

has integrable (square root) singularities at x - a and x = b. In solving 

the integral equations numerically it is assumed that 

ilf(t) = G(t)[(t-a)(b-t)]‘ 1/2 , a<t<b , (63) 

where 6 is bounded in a<t<b. Thus the index of (62) is +1 and the solu- 
tion is subject to the following additional condition 


b 

r (a,0) = v = - *(t)dt . 
u J a 


The integrals involving ^(t) may be evaluated by using the Gauss-Chebyshev 
integration formulas given in [9] after normalizing the interval (a,b) as 

S = (2x-a-b)/(b-a) , t = (2t-a-b)/(b-a) . (65) 


To evaluate the integrals involving $(r) first the variables y and r 
are changed as 
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y = h - n , 


r = h - p 


( 66 ) 


and then the definition of <f> and the kernels are extended into -h<n<0, 
-h<p<0. By defining now 

1/2 

<j»(r) = 4>(h-p) = <(> 0 (p) = g(p)[h 2 -p 2 ] , (67) 

the Gauss-Chebyshev integration formulas may then be used to evaluate the 
integrals [9,8]. 

After the density functions are obtained the stress intensity factor 
k(b) at the crack tip and the strength of the compressive stress singu- 
larity k (a ) at x = a may be obtained as 


k(b) = lim /2(b-x) a (x,0) 
x-+b yy 

= " i+7 lim /2 ( b " x ' *( x ) » 

' K x-*b 


( 68 ) 


k (a ) = lim /2(x-a) a (x,0) 
x-*a yy 


= lim /2(x-a) *(x) . 

X ^3 


(69) 


Figure 2 shows some sample results for the crack surface displacement 
which is obtained from 


v(x,0) = - 


\p( t ) dt , a<x<b , 


(70) 


The contact stress distribution for various values b/h and a/h is shown in 
Figure 3. In this figure the stress a yy (x,0) is normalized with respect to 
the average contact stress a Q given by 

r* 

o 


1 r f v (Xi0)dx • 


(71) 
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The stress a yy( x »0) is obtained from (10) which is valid for 0<x<®. The 
corresponding stress intensity factors are shown in Table 1 which also 
shows the resultant compressive wedging force per unit thickness calculated 
from 


P 


y 


/-a 

^ ° yy (x,0)dx = ao 0 


(72) 


TABLE 1 

Stress Intensity Factor for a Flat Wedge: v Q /h = 0.1 


b/h 

a/h 

V 102 

k(b) 

P y -b/h 377 

1.0 

0.1 

1.60 

5.928 

1.0 

0.2 

1.86 

5.983 

1.0 

0.3 

2.13 

6.051 

1.0 

0.37 

2.35 

6.149 

2.0 

0.1 

0.494 

4.670 

3.0 

0.1 

0.209 

4.272 

4.0 

0.1 

0.107 

4.055 


From Figure 3 it may be observed that for a fixed wedge thickness 2v Q 
and fixed crack length b, as penetration distance a increases the contact 
stress Oyy(x,0) at x = 0 decreases, and at a certain value of a = a! it be- 
comes zero. For a>a-j the procedure described above to solve the problem 
is not valid and the problem must be formulated in terms of a system of 
three integral equations as indicated by (17-19). In the example shown 
by Figure 3, v o = 0.1h, b = h, and the value a-| , at which the wedge starts 
separating from the strip is found to be a-j = 0.37b. 

As a second example we consider loading through a rigid frictionless 
wedge with a constant slope (see the insert in Figure 4). In this case 
the crack surface tractions ay y (x,0), a<x<b, and o xx (0,y), 0<y<h, are again 
zero and the integral equations become (see equations (9) and (10)) 
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where 



•h 

o 


H n (y»r)<|>(r)dr + 
H 21 (x,r)<j>(r)dr + 


■b 

H 12 (y.t)*(t)dt = 0 , 0<y<h , (73) 


b ra 

H 22 (x,t)\j;(t)dt = - H 22 (x,t)f(t)dt , 

^ ■ A 


a<x<b , (74) 


f(x) = v(x,+0) = -f 0 , 0<x<a . (75) 

Because of the smooth contact at x = a, the density function ij>(x) is of the 
following form 


*(x) = -f Q + G(x)[(x-a)/(b-x)] 1/2 . (76) 

In solving the integral equations numerically we also modify the definition 
of $(r) as 


4>(r) = $(h-p) = $ q (p) = f o + g(p)(h 2 -p 2 )^ 2 


(77) 


Note thatwitk (76) and (77) the condition (58) is automatically satisfied. 

In this problem the known loading parameter is either the wedge driv- 
ing force P Q or the wedge thickness 6 at the strip end, i.e., 2v(0,0). The 
wedge thickness is calculated from 


5 = 2v(0,n) = 2(af 0 


r b 

'!>( t)dt) 
J a 


(78) 


Equations (71) and (72) are still valid giving the average contact stress 
o Q and the wedging force P y , respectively. Once P^ is found, the wedge 
driving force may be obtained from 


P 0 = 2P y sin(tan"f 0 ) = 2P y f 0 


(79) 
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In this problem the contact distance a is unknown and is determined either 
from the equilibrium condition (72) or from (78). However, the prob 1 em is 
more conveniently solved by specifying a, solving the integral equations, 
and then determining the corresponding load P Q or the thickness 6. Note 
that in this smooth contact problem the index of the integral equation (10) 
is zero and therefore no additional condition is required for a unique 
sol ution. 

For a fixed crack length b = h and wedge slope, f Q = 0.01, Figure 4 
shows the calculated values of the stress intensity factor as a function 
of 6, the crack opening at x = 0. Figure 5 shows a sample result for the 
distribution of the contact stress, a which is again normalized with re- 
spect to the average stress o 0 given by (71). Figure 6 shows the crack 
surface displacement for various values of the crack length. Note that 
in order to have the same contact distance a = 0.1h, in a long crack the 
penetration of the wedge (or the crack opening at x = 0) needs to be 
much greater than that needed for a short crack. 

The calculated stress intensity factors for a rigid wedge with a 
constant slope f are summarized in Table 2. As indicated before, in the 
numerical analysis the contact distance a is specified and then the corre- 
sponding values 6 = 2v(0,0) and the loads (per unit thickness) P y and P Q 
are calculated. The stress intensity factor k(b) defined by (68) is shown 
in the last column. The particular normalization factor is introduced in 
order to compare the results with that of Ref. [2] and the beam theory. 

This comparison is shown in Figure 7. Using the classical beam theory, 
simple energy balance consideration would give the stress intensity factor 
as follows: 

k(b) = A2 P y b/h 3/2 (80) 

where P y is a concentrated wedging force (per unit thickness) applied at 
the end x = 0. If the strain energy due to shear deformations is also taken 
into account, (80) is slightly modified and becomes 

1 /2 

k(b) • m (p y b/h 3/2 )(i + (si) 
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TABLE 2 

Stress Intensity Factor for a Rigid Wedge with Constant Slope 


b/h 

a/h 

f 

0 

:x|ch 

o 

ro 

o+<)p„ 
y i n 3 

OMP n 

0 in s 

. Mb) 

4yh 10 

4^h 10 

Pyb/h 3/2 

1.0 

0.05 

0.01 

0.597 

1.31 

2.62 

5.918 

1.0 

0.1 

0.01 

1.078 

2.53 

5.05 

5.926 

1.0 

0.15 

0.01 

1.458 

3.57 

7.14 

5.943 

1.0 

0.2 

0.01 

1.735 

4.38 

8.77 

5.950 

1.0 

0.4 

0.01 

2.166 

5.84 

11.7 

6.065 

1.0 

0.1 

0.04 

4.314 

1.01 

80.8 

5.926 

2.0 

0.1 

0.01 

2.838 

2.31 

4.62 

4.597 

3.0 

0.1 

0.01 

5.791 

2.12 

4.23 

4.150 

4.0 

0.1 

0.01 

9.962 

1.94 

3.38 

3.912 

5.0 

0.1 

0.01 

15.33 

1.79 

3.53 

3.715 


Figure 7 shows the results obtained from the wedge solution (with 
a = 0.1h), Ref. [2], (80), and (81). The results given in [2] were calcu- 
lated by using a boundary collocation method for a strip which has a fin- 
ite length. Considering the differences in loading and geometry the 
agreement is reasonable. 

Finally, a sample result for the crack opening displacement in a 
strip subjected to a distributed wedge force is shown in Figure 8. Here 
the crack surface traction is assumed to be 

o yy (x,+0) = - /d 2 -(x-0) 2 . (82) 

This is an elliptic distribution which may be representative of smooth 
contact problems with elastic wedges. In the example the constants are 
assumed to be 

b s 1 .143h , D s 0.524h , d * 0.2h . (83) 
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The b/h ratio considered here corresponds to the standard compact tension 
specimen. For this example the stress intensity factor is found as 


= 5.971 

P y b/h J '^ 


„ _ 2 h D 
p o ' (T cP P y 


(84) 


where is the resultant force per unit thickness. 

Using the analysis given in this paper one may easily calculate a 
"Green's function" for the stress intensity factor k ( b ) which can be used 
to solve the edge crack problem under an arbitrarily distributed crack 
surface pressure. To do this in (9) and (10) we substitute c = 0, d s h, 
cr xx = (0,y) and assume that 


(-1 for x = x< , 

a vv (x,0) = < (85) 

yy [ 0 for x = x-, j = l,..,i-l, i+l,..,n 

J 

where x./b, i = l,..,n are the appropriate collocation points in (0,1) used 
in the solution of the singular integral equations. Solving the problem 
for the loading (85) the stress intensity factor is obtained as a function 
of x.j in the following form: 


k(b,x i ) = G(x i ) , i = 1 , . . r n . (86) 

If the crack surface pressure is now specified as 

ayy(x,0) = -f(x) , 0<x<b , (87) 

then the stress intensity factor may be obtained as follows: 

k(b) = l G(x.)f(x.) . (88) 

i*l 1 1 

This procedure follows from the method of solution used in solving the 
singular integral equations. It may be noted that G(x^ ) also corresponds 
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to the stress intensity factor due to a concentrated wedging force of 
unit magnitude applied at x = xi. Table 3 shows an example obtained for 
b = h. 


TABLE 3 

Green's Function for a Symmetric Crack 
in a Semi-infinite Strip, b = h 


i 

1 

2 

3 

4 

5 

xi/b 


0.013 

0.036 

0.070 

0.114 

GUi)/(b/h3/2) 

5.915 

5.919 

5.923 

5.930 

5.951 


6 

7 

8 

9 

10 

11 

12 

0.167 

0.228 

0.295 

0.368- 


0.519 

0.595 

5.986 

6.063 

6.169 

6.329 

6.495 

6.690 

7.009 



14 

15 

16 

17 

18 

SHI 

20 




mm 

! 1 

0.948 


0.994 




ESS 

ESS 

15.75 


44.95 
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Figure 1 Geometry of the 


cracked strip 






gure 2 Crack surface displacement in a strip loaded by 
flat wedge, v 0 = O.lh 
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Figure 3 Distribution of the contact stress in a strip loaded 
by a flat wedge, b=h, v o = 0.1h 
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Stress intensity factor in a strip loaded by a rigid 
wedge with constant slope. av(x,+C)/3x = -f 0 = -0 .01 , 
b = h, S = 2v(0,0) (crack opening at x = 0) 
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Figure 7 


b/h 

Stress intensity factor in a strip loaded by a wedge 
with constant slope, f o = 0.01, a = 0.1h 
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